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Abstract 

We study a conjecture by Fendley, Ludwig and Saleur for the nonlinear conduc- 
tance in the boundary sine-Gordon model. They have calculated the perturba- 
tive series of twisted partition functions, which require particular (unphysical) 
imaginary values of the bias, by applying the tools of Jack symmetric functions 
to the "log-sine" Coulomb gas on a circle. We have analyzed the conjectured 
relation between the analytically continued free energy and the nonlinear con- 
ductance in various limits. We confirm the conjecture for weak and strong 
tunneling, in the classical regime, and in the zero temperature limit. We also 
shed light on this special variant of the Imi^-method and compare it with the 
real-time Keldysh approach. In addition, wc address the issue of quantum sta- 
tistical fluctuations. 
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1. Introduction 

Quantum Brownian motion (QBM) of a particle or charge in a tilted peri- 
odic potential is one of the basic subjects of nonequilibrium quantum-statistical 



i: 



condensed matter [1| . Of particular interest is a duality symmetry between the 



mechanics and a key_ model for a multitude of quantum transport phenomena in 

symrr 

weak- and strong-binding representations of the model 0, Q . There is a close 
correspondence of the QBM model with seemingly different models: (i) a resis- 
tively shunted Josephson junction jij], (ii) a quantum impurities in a ID wire 
(iii) tunneling of edge currents through constrictions in fractional quantum-Hall 
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devices |6(, and (iv) a coherent one-channel conductor in a resistive electro- 
magnetic environment 0j- In fact, all these models are directly related to each 
other. 

The calculation of the dynamics may be done with the use of the Feynman- 
Vernon influence functional method Q or the equivalent nonequilibrium Keldysh 
technique. In this theory, the reduced density matrix (RDM) is expressed as 
a double path sum over all possible paths on the RDM-plane. The path sum 
problem is formally equivalent to the grand-canonical sum of a charge gas with 
complex-valued interactions. Of particular interest is the scaling or field theory 
limit. Until now the respective real-time Coulomb gas has not been solved 
analytically barring special cases. 

One may also consider the system in imaginary time compactificd on a circle 
of length 1/T, where T is the temperature. In the field theory limit, the charges 
of the respective Coulomb gas interact with a "log-sine" long range interaction. 
The compactified Coulomb gas has been solved analytically for particular un- 
physical values of the bias with the tool of Jack symmetric functions The 
consideration of nonequilibrium quantum transport requires (i) analytic con- 
tinuation of the calculated partition function to the physical bias, and (ii) a 
conjectured relation between the analytically continued partition function and 
transport quantities Here we analyze this route in different parameter 

regimes, and we compare the findings with results available by different meth- 
ods. 

Our paper is organized as follows. In Section [2] we introduce several related 
models, which may describe nonequilibrium quantum transport of a particle or 
a charge, Josephson junction dynamics, scattering off an impurity embedded 
in a Luttinger liquid, and conductance of one-channel conductors in a resistive 
electromagnetic environment. In Section [3] we review the solution of the "log- 
sine" Coulomb gas on a circle with winding numbers. In Section [4] we study the 
conjecture for the nonlinear conductance in the weak-tunneling regime, in the 
classical regime, and in the zero temperature limit. In Section [5] we compare 
the results of the imaginary time approach with those of the real-time Keldysh 
method. We also address the issue of statistical fluctuations. Section [5] is 
devoted to concluding remarks. 

2. Generic transport model 

The boundary sine-Gordon (BSC) model is a one-dimensional quantum field 
with gapless bulk excitations and a sinusoidal boundary interaction. The model 
is integrable and can be treated by a variety of powerful techniques. It con- 
stitutes a generic quantum transport model which displays a wide range of 
interesting characteristics, and it is a realistic model to describe seemingly dif- 
ferent physical situations. In fact, all these models are closely related to each 
other. They have in common that a nonlinear quantum degree of freedom is 
coupled to a harmonic field. 
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2.1. Nonlinear quantum Brownian motion 

A quantum Brownian particle moving in a periodic potential is a key model 
for many transport phenomena in condensed matter. Of particular interest is 
the case of strict Ohmic dissipation in which the regions of weak and strong 
corrugation are related by an exact self-duality. 

2.2. The dual representations 

The dissipativc dynamics of a particle moving in a tilted, sinusoidally corru- 
gated potential of periodicity length X may be described by the translational- 
invariant Hamiltonian with bilinear environmental coupling^] 



P 2 / 2ttX\ X 

t ^ 



Xq J V w ra x r 



(i) 



The "weak-binding" form (|T|) is the convenient representation for the pertur- 
bative expansion of transport quantities in the corrugation strength U . Of 
particular interest is the normalized mobility M. = n/fio (particle transport) or 
normalized conductance Q = G/Gq (charge transport), where no and Go are 
mobility and conductance in the absence of the corrugation. We have Ai = Q. 
In the sequel we use the language of charge transport. 

In the tight-binding regime T/lj , e TB /uj -C 1 <C Uq/lu , where w 
is the distance of low-lying levels in the individual well regimes, and e TB is 
the potential drop between neighboring wells, only the lowest energy band is 
occupied. Then the system is effectively reduced to a tight-binding lattice. The 
tight-binding or tunneling Hamiltonian with transfer matrix element A and 
translational invariant coupling is 



#tb = --^^2(\n)(n+l\ + \n+l)(n\) - e TB Y 



y 



(2) 



where Y Q is the TB lattice constant, and 



Y = Y J2n\n)(n\ (3) 

n 

is the position in TB representation. Since wc postulate the same tilting force 
in both models, the potential drops between neighboring sites are related by 



-wb 



/X = e TB /Y . (4) 



Throughout this paper we put h = fc B = 1. 
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In the scaling or field theory limit, the spectral densities of the environmental 
couplings, 

(5J 

G TB( W ) = EaAia<5(w-^ Q ) , 

are of strict Ohmic form 

Gwb(^) = 2A'wbw , 

(6) 

Gtb(w) = 2Ktb <*> ■ 

The dimensionless so-called Kondo parameters Awb and Atb are related to 
the Ohmic viscosity 77 by 

In the Ohmic regime ([6]), the TB and WB model ([I]) and ([2]) are related by an 
exact self-dual symmetry, as noted first by Schmid for the equilibrium density 
matrix at T — 0- Later, the duality was proven to hold also at finite T and 
for real time dynamics Q . In the self-dual mapping we have 

YJX = Atb ee K = 1/A WB , (8) 

and the normalized conductance in the WB model ([T]) is directly related to the 
corresponding quantity in the associated TB model ([2]), 

G WB (T, e, K) = 1 - G TB (T, e/K, l/K) . (9) 

The self-duality relation ^ exchanges weak and strong corrugation, as well 
as the regimes of weak and strong environmental coupling. The explicit calcu- 
lation of ([9]) may be executed with the use of the Feynman- Vernon influence 
functional method |8[, or directly by unitary transformations of the WB and 
TB Hamiltonians 

Computations are conveniently performed in the discrete TB model. With 
use of the self-duality, the results can then easily be transferred to the dual WB 
model (see below). 

With regard to the partition function, we consider the system in imaginary 
time compactified on a circle of length 1/T. The partion functions -Ztb and 
ZyjB can be expanded in even powers of A and Uo, respectively. The pertur- 
bative series of the partition function is equivalent to the grand-canonical sum 
of a one-dimensional gas of positive and negative charges £ = ±1 with overall 
neutrality. The charges represent the forward and backward moves in the TB 
or WB potential. All effects of the environmental coupling are described by the 
imaginary-time charge-charge correlator 

/7^ e i27rg(r)/g e -i27rg(0)/ 9o \ _ g - (27r/ 9o ) 2 ( [ ?(0)-«(r) ] g(0) > 

(10) 
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where (• • • ) means quantum statistical equilibrium average. The spectral rep- 
resentation of the charge interaction is 



/>oo 

W{t) = / dw 
Jo 



G(ijj) cosh(^) - cosh[u;(i - r))] 



sinh( 



2T > 



In the scaling limit (J6j> we have 

W{t) = 2K In 



ttT 



sin(7rTT) 



(11) 



(12) 



where ui c is the UV- frequency cutoff in the spectral density G(ui). 

While in a two-state system the charges alternate, in our case of infinitely 
many states the charges are unordered. In the term of order A 2n or [/q™, we have 
n positive and n negative charges, and we must integrate over the locations of 
the charges in the interval 1/T. The grand-canonical series of the TB partition 
function is 



Z 



TB 



(A/2) 



2n rl/T n 



JO i = 1 Jo . ! 



xexpJ E [W{r m -T k ) + W{r' m -r' k )] 



^ m> k — 1 



(13) 



exp<^- E W(\r m -r' k \) 

[ k,m=l 

x exp|-eE( T fc- r fc)l ' 



k k = l ) 

and a similar expression holds for 2\vb- 

3. Josephson junction dynamics 
According to the dc Josephson effect, a super current of Cooper pairs flows 
through a superconducting tunnel junction. At finite voltage V x , there is a 
supercurrent only if the excess energy 2eV x can drain off to the electromagnetic 
environment. Quasiparticle excitations can be disregarded when the voltage V x 
is small compared to the gap voltage. For weak Josephson coupling, E 3 <C E c = 
2e 2 /C, the charge on the junction is fairly well defined. The relevant minimal 
model of the junction-plus-environment entity is 12 1 



H 



2C 

N 

E 

a=l 



Ej COS 1p 



a 2 
2C\ 



1 1 

~c^2L r 



(14) 



The first term represents the charging energy, the second term accounts for 
tunneling of Cooper pairs through the junction, and the third term describes 



5 



the coupling to the electromagnetic environment, which is modeled by a set of 
LC-circuits, and it also includes the applied voltage. In the limit N — > oo, the 
capacitances {C a } and inductances {L a } provide an arbitrary linear admittance 



N 



Y(w) = Z~\u}) = / dte~ iwt V 

J a = l 



cos(cj Q t) 

in 



where uj a = 1/ \/L a C a . The total impedance seen by the junction is 



\ujC + Y(uj) 



(15) 



(16) 



After tracing out the electromagnetic modes, one finds for the imaginary-time 
phase correlator in thermal equilibrium 



where W^(r) has the spectral representation (| 1 1 1) with 

G^{u) = 2 cj , 



(0)> = e -Wv,(r) 



(17) 



(18) 



and where i?q = = ir/2e 2 is the resistance quantum for Cooper pairs. 

For a resistive electromagnetic environment, ReZ t (o;)/i?Q = p/[l + (w/w R ) 2 ] 
with p = R/Rq, and lo r = 1/[Z(Q)C] = E c /(irp), we get 



W^t) = 2pln sin^IY) 
L 7TJ 



2pC 



(19) 



The first term in Eq. (|19[) describes the resistive environment in the scaling limit 
[ cf. Eq. (jT2j ] , and the quantity £ encapsulates details of the total impedance 
at high frequencies 13( , 



c 



*( 1+ &)-*w- h (&) 

ReZ t (w) 1 



'dw 



pRq 



1 + (w/w R ) 2 



(20) 



The second term in (flU]) leads to adiabatic renormalization of the tunneling 
coupling 

Ej -> = £je-' £ . (21) 

The model (IT?|) with (jT9")l directly corresponds to the model © with ©. 
The correspondence relations are 



2eK 



K = p 



(22) 



In the opposite limit of large Josephson coupling, Ej 3> E c , the phase of 
the BCS condensate is fairly well defined and localized in one of the wells of 



G 



the Josephson potential. Then it is pertinent to turn to the Bloch band rep- 
resentation of the Josephson junction and to choose a basis of discrete phase 
states, 

i> = 2n^n\n)(n\ , e^ 2e = £>><n+l| , (23) 

n n 

where Q is the quasi-chargc. Assuming that the junction dynamics is confined 
to the lowest energy band of width U , we arrive at the Hamiltonian (3. [l2| 



H 



U cos 



2ttQ 
~2e~ 



2n 



2C„ 



1 1 



1> T 



(24) 



For a current-biased junction we have V K = RqI x . Since charge and phase 
arc interchanged when turning from the charge representation (|14p to the phase 
representation (|24l) . the environmental effects are now captured by the charge 
autocorrelation function Wq(t) = (ir/e) 2 ([Q(0) - Q(t)]Q(0)). The role of 
Zt(w)/i?Q in (|18p is now taken by RqY(lu), where Y(us) is the admittance. 
Thus we obtain for a strictly resistive environment 



T%(r) = -ln 



-| sin(^Tr) 



(25) 



where oj r is the required cutoff. 



2. 4- Quantum impurity in a Luttinger liquid 

The ground state of ID interacting spinlcss electrons is a Tomonaga-Luttingcr 
liquid (TLL), which is characterized by a gapless collective sound mode. The 
low-energy modes of the ID interacting electron liquid is conveniently treated in 
the framework of bosonization. An extremely useful field-theoretical formula- 
tion of a TLL has been given byHaldane 0. The creation operator for spinless 
fermions may be expressed in terms of bosonic phase field operators <j)(x, t) and 
9(x,t) as %p'(x) cx X)nodd ex P{i n [' c F a:; + V^^i x )}} °xp[i v / 7T</)(a;)]. At long wave 
lengths only the terms n = ±1 are important. They represent the right- and 
left-moving parts of the electron field. The bosonic field operators obey the 
equal-time commutation relation [(j)(x,t), 6(x',t)] = — isgn(a; — x')/2. 

The generic interaction of the TLL universality class is effectively given by a 
(S-potcntial. This yields the interaction term H 1 = U J dxdy p(x) 5(x — y) p(y). 
In the TLL, the electron interaction is characterized by a single dimensionless 
parameter g = l/-v/I + U /ttv F: where v F is the Fermi velocity. Excitations of 
the liquid are described by the generic harmonic Hamiltonian 

H L (g) = \\^ [{dxOf/g + g{d x 4>) 2 ] , (26) 
where v = v F / g is the sound velocity. 
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Consider now a single strong point-like impurity, or equivalently a weak link 
or tunnel junction. A hopping or tunneling term which transfers electrons across 
the impurity or weak link may be included in the original electronic Hamiltonian 
by adding H[ = — A[ip^(x = + ) tp(x = 0~) + h. c. ]. Expressed in terms of the 
bosonic fields, the hopping term induces a jump of the (/>-field at the impurity, 
4> = ^[4>{ x = + ) — <t>{x = - )]. In the (^-representation, the impurity is then 
described by 

H[(~4>) = - Acos[20F0 + eV a t] . (27) 

Here we have included a voltage term. The model is rounded off with addition of 
the harmonic liquid in the right/left (±) lead. We then arrive at the weak-link 
Hamiltonian (</>-model) [H[ 

= H Li+ (g) + H L _(g) + H{@) . (28) 

Since the nonlinear degree of freedom is only at x = 0, we may inte- 
grate out the fluctuations of 4>(x) away from x = 0. With the Fourier ansatz 
0(.t, t) = TJ2 n 4>(x, v n ) e~ [1/ ™ T , where v n = 2irTn is a Matsubara frequency, the 
Euclidean action is minimized when (j)(x,v n ) = <p(v n ) exp(— \v n x\/v). The re- 
sulting Euclidean influence action, which includes all effects of the excitations in 
the leads on the tunneling degree of freedom, is S^((f>) = gT^2 n \v n \ \4>(v n )\ 2 . 
Average with the weight function e _s i=fi® yields for the charge correlator 

7 - e i2^(r) e -i2V5F0(O)\ = e -W-(r) (2g) 

with the charge interaction 



= - ln[^ S in(7rTr)' . 
A single weak impurity is modeled by the Hamiltonian 



(30) 



H sc = / dxU(x)^(x)^(x) , (31) 



where U(x) is the scattering potential. For a point-like scattcrcr at x — 0, 
the 2fc F -backscattering contribution takes the form H sc (9) = — U cos[2y / 7r 6], 
where 6 = 0(0). With the lead and voltage term added, the weak-impurity or 
strong-tunneling Hamiltonian (0-model) is 

H e = -U cos[2^6] + eV6/^ + H L (g) . (32) 

One may again integrate out the Luttingcr modes away from the scatterer. This 
leads to the influence action for the scattering mode S^(0) = {^/g)T^2 n \v n \ \8(v n )\ 2 , 
and the charge interaction in the Coulomb gas representation takes the form 

W„(t) = 2g ln[^ sm^Tr)] . (33) 
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The (/>-model is self-dual to the 0-model, and vice versa. In the mapping the 
regimes of a weak and a strong barrier are interchanged. We also see that the 
4>- and ^-representations of the impurity model directly correspond to the above 
TB- and WB-representations of the Brownian particle in the washboard poten- 
tial. The Kondo parameter K is related to the electronic interaction parameter 
g by 

g = 1/K . (34) 

Since the correspondences are valid for the Hamiltonians and the effective ac- 
tions [H, they hold not only for transport quantities like the current but also 
for all quantum statistical fluctuations. 

In addtion, the correspondences also hold for the full counting statistics of 
tunneling of edge currents in fractional quantum Hall (FQH) systems, where 
the fractional filling factor v takes the role of g (see e.g. Ref. [6(). 

2.5. One-channel coherent conductor in a resistive electromagnetic environment 
The ID fermion-boson correspondence suggests that a coherent one-channel 
conductor under influence of a resistive electromagnetic environment can be 
mapped on a quantum impurity embedded in a TLL. The mapping has been 
analyzed by Safi and Saleur Q. 

A mesoscopic conductor in contact with an electromagnetic environment 
forms a quantum system violating Ohm's law. Transfer of energy from the elec- 
trons to the environment leads to dynamical Coulomb-blockade, which reduces 
the current at low voltage. The Luttingcr parameter for a coherent conductor in 
the absence of the electromagnetic environment is g = I . 

Following the treatment given in the preceding subsection, and augmenting 
the phase in the tunneling term (|27j) with the phase ip of the electromagnetic 
environment, the relevant weak-tunneling Hamiltonian is found to read 

H wt = - Acos[2V^0 - <p + eVJ] 

(35) 

+ H h+ (g = 1) + H L _(g = 1 ) + H env [Q,ip] . 

It is convenient to combine the electronic phase with the electromagnetic phase 
into the auxiliary phase \ = <j> — ip. Following the route sketched in the 
previous subsection, one finds for the phase correlation function for a resistive 
environment ( [x(0) — x( T )] x(0) ) = W / y( t ) w ^ tn 

Wrfr) = 2(1 + p) In [ ^ sin(TrTT) ] + 2pC ■ (36) 

Here p = R/Rk, and i?x = 27r/e 2 is the resistance quantum. 

The Hamiltonian of the corresponding strong-tunneling conductor is 

H st = - U cos[2^e} + H L (g = 1) 

-(eV a + p)0/V^+ H env [Q,cp] , 

with V a the applied voltage and ip the fluctuating phase of the environment. 
Upon integrating out the electronic modes away from x = 0, and also the 
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environmental modes for a strictly resistive environment, the effective action 
takes the form S^{8) = (1 + p)T \v n \ \9{v n )\ 2 . Evidently, if we turn to the 
Coulomb gas representation, this leads to the charge interaction 

W ¥ (r) = jl- In [ ^ sin(^Tr) 1 + 2(/p . (38) 

Again, we have added both in ([31))) and in (J3HJ) a term which captures the 
details of the resistive environment at high frequencies as given in Eq. (|20[) . This 
term leads to an adiabatic renormalization of the couplings A and Uq, as in Eq. 

Thus we find that the two representations of the coherent conductor are self- 
dual to each other. Furthermore, they can be mapped on the models described 
before. In the mapping to the Brownian particle model and to the quantum 
impurity system we have 

K = 1/g = 1 + a . (39) 

These correspondence relations hold both for the weak- and strong-tunneling 
representations. 



3. Coulomb gas representation 

In the sequel, we choose from the above diverse models the tight-binding 
or weak-tunneling model ©■ With the duality mapping flSJ) and © and the 
above correspondences between the various models, subsequent findings can be 
transferred easily to the weak-binding regime and to both regimes of the other 
models. 

The change of integration variables, itj = maps the Coulomb gas (|13p 

with (|12p on the unit circle. The pcrturbativc scries of the partition function 
(TT51) may be written as a power series in the effective fugacity 



A /2ttT 
X ~ 2T 

The tight-binding series takes the form 



0J C 



K 



(40) 



Z TB (x,p) = 1+ V x 2n l 2n (p) , 



71=1 



(41) 



where 



-2Kn 



(n!) 2 

Ili<j sin(^) sin(^P) 



n 



i= 1 

f Ua—U. 



du 4 du'i 
~2ir ~2ir 



2 A' 



(42) 
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and where p represents the scaled physical bias, 



p = iq with q = -^—f ■ ( 43 ) 



Changing integration variables Zj = c 1Ui and z[ = e 1Ui we get 



[A(z)A{z)] K [A{z')A{z')] K (zx~-Zn 

n <lfc [(i-M*)(i-«] 



(44) 



where A(z) = Oi<fc( z i ~~ z k) IS the n- variable Vandermonde determinant. 

Unfortunately, the multiple integrals in (|42|) or (|44|) can not be evaluated 
for general if in the regime < K < i and for general complex p. To advance, 
we must put K rational and p integer. Integer p may be regarded as a winding 
number due to a magnetic charge located at the origin. For rational K , we may 
expand the integrand in terms of Jack polynomials 0, [TBI ] 



(l - nsj) 



bx(K)P x (r, K) P x (s, K) . (45) 



K 



Here, the function P\(r,K) is a symmetric polynomial in the set of variables 
{ti, t%,- ■ ■ , r„), and A = (Ax, A2, ■ • • , A„) is a partition of an integer, Ai < 
A 2 < • • • < A n . For positive integer p, we have 

( Zl ---z n yP x (z,K) = P x+p (z,K), (46) 

where A +p means the partition A where p columns of length n have been added 
in the respective Young tableau. The multiple integrals in (j44]l can be executed 
by use of the orthogonality relation of the Jack polynomials, since it involves 
the Vandermonde determinant. In the end, one arrives at the totally ordered 
n-fold series expression 

00 rn n ni2 n 

x 2n ( P )= J2 E ■■■ E II^K) (47) 

m n — rrin — l— mi— j — 1 



with 



&j (to) 



1 T{jK + to) 



T 2 (K) T(l - K+jK + m) 
T(jK + p + to) 



(48) 



T(l - K + jK +p + m) 

It has been suggested in Ref. [h| to use (gT|) with g7]) and (gSJ) to define 
Z TB (x,p) for complex p. Furthermore, it has been conjectured that this is 
the unique analytic continuation. Below, we confirm this conjecture in various 
limiting cases where analytic expressions are also available by different methods. 
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For instance, we find from the significant limit \p\ — > oo that the free energgy 
F(x,p) = —T\nZ(x,p) has the proper behavior in all orders of A, as T — > 0. 
Consider now the cumulant expansion 

oo 

\nZ TB (x,p) = ]Tx 2n C 2n (p) . (49) 

n = l 

The first three cumulant coefficients are 
C 2 (p) = I 2 (p) , 

dip) = li(p) - i2f(p) , (50) 

d(p) = Mp)-Mp)Mp) + \iI{p) ■ 

For the analysis of the limit \p\ — > oo in the cumulant expressions, it will be 
pivotal to rearrange the multiple sums in C2n (p) into totally ordered sums as in 
(|47]) . This leads to the split-up 

n 

?("»)/ 

m — 1 

i(m). 



C an (p) = J2 C ^\p) . (51) 

m — 1 

where (p) is an m-fold ordered sum. The first cumulant is 



C 2 (P) = E e i(j) : ( 52 ) 

3 = 1 



The second cumulant has two contributions, 
C1 X) (P) = 5 f>?(j), 

3 = 
oo 3 

= EE [ea(i)-ei(j)]ei(AO 

j = k = 

The three contributions of the third cumulant are 

ci 1] {p) = \ £e?(j), 

6 3=0 

OO 

7(2)/ 



(53) 



j>k = o (54) 



4 3) (p) = E {[e 3 (i)-ei(i)WfcKM 

f = o 

-2[e 2 (i)-e 1 (i)]e 1 (fc)e 1 w} 



j > fc > £ = o 



The integrals in the expression (J42J) , which defines the perturbative coeffi- 
cients l2n(p), diverge at short distances when K > 1/2. Correspondingly, the 
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multiple series ([4"Tf for the coefficient l2n(p) diverges in the regime K > 1/2 for 
all n. Interestingly, the coefficient C2n{p) in the cumulant series (|49| is regular 
for A"- values even when T 2n {p) is singular. This is because some of the diver- 
gences in the Coulomb integrals are cancelled when taking the connected part 
Cin (p) ■ The analysis of the ordered m-fold sums (p) in dHJ) shows that the 
cumulant coefficient C 2n (p) is nonsingular in the range < K < 1 — (see 
below for details). 



4. Nonlinear conductance 

In Ref. (lol | a conjecture was made which relates the nonlinear conductance 
directly to the partition function Z(x,p) or to the free energy ^(x^p). The 
conjecture for the conductance of the TB model @ is 



r , /rr s K d ( Z TB (x,-p) 
Gtb [x, e/T) = — x— In 



2p dx \ Ztb{x,p) 
K x d 



2pTdx 



(Ftb{x,p) -T T b(x, -p)) , 



(55) 



where p is given in Eq. (|43|) . The conjecture is based on the fact that the free 
energy is real only for p integer, and that the continuation of C2n{p) is not even 
in p, so that the free energy aquires an imaginary part when p is complex. 
The respective relation for the WB model (H} is 

Swb(x, e/T) = 1 - -L d ln ( ^-P) \ (56) 
2Kp dx \ Z WB {x,p) J 

where p = p/K = ie/(2iTKT). 

Use of the perturbative cumulant expansion (|49)) in the relation (|55|) yields 
the weak-tunneling series of the conductance, 

oo 

g T B(x,e/T) = J2G n {x,e/T). (57) 

n = l 

The perturbative tunneling contribution of order n is 

9 n (x,V/T) = ^ n [C 2n (-iq) - C 2n (iq)]x 2n , (58) 
l q 1 J 

where q = e/ (2ttT). We now prove the conjecture in diverse parameter regimes. 

4-1- Strong barrier limit 

The series that defines C2(p) can be summed in analytic form 



r( , 1 ^ T(K + j)T(K + p + j) 

2[P) T*(K) f^ Q r(l+j)T(l+p+j) 

sm(KTT)T(l-2K) T(K + p) 

7T r(l- K+p) ' 



(59) 
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With this form the perturbative conductance (|58p in order A 2 is found as 

Qx = K /AN 2 (2nT\ 2K sinhp |F(A + i^)| 2 (6Q) 



2TJ V / # r(2A) 

Correspondingly, the perturbative conductance through a voltage-biased 
Josephson junction arising from incoherent Cooper pair tunneling is 

E' L \ 2 (2^\ 2p sinh(f ) |r(p + i^)| 2 (g 



2Tj V K J r(2p) 

The expressions (|60[) and ()Q 1 [) are exactly the same as the ones obtained from 
the non-equilibrium real-time approach 0, [l3[ . 

Let us study the equivalence in some detail. In the weak tunneling limit, 
the normalized conductance is given in terms of the golden rule transition rates 
fcj 1 representing incoherent forward/backward tunneling transitions to adja- 
cent wells (see Section [5] for notation). We have 

Oi = § (*# - *#) , (62) 

where q is given in (|4"5| . The rates may be calculated directly from the standard 
real-time double-path representation in order A 2 , which is pictorially sketched 
in the right diagrams of Fig. HQ Alternatively, one may calculate the rates with 
the standard Im ^-method Here, the free energy of order A 2 , Ti(=fp), 
(see left diagrams of Fig. [1]) is analytically continued to a complex value. The 
imaginary part of the free energy is obtained by distorting the integration path 
of the charge correlation integral from imaginary to real time. In this way we 
get 

«=<f> - . — d» — w . (63) 

Taking either road we obtain the unique expression [l[ 

= x 2 T [ e- 1 K * Ji (K, tp)+c 1 K * J\ (K, ± p) ] , (64) 
where J± (A, p) is the real-time noise correlation integral 



1 f°° 

Ji(K,p) = - du 



(2sinhu) 2 ^ 
T{1-2K)T{K + p) 
2irT{l -K+p) 

Upon combining (|52"j) with we get 

2A 



(65) 



— x 2 sin(A7r) [ J x (A, -p) - J x (A, p) ] , (66) 
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Figure 1: Imaginary-time charge diagrams representing Z\{^-p) or the connected part T\(py-V) 
(left), and analytic continuation to real-time forward/backward path (right). Diagram (a) 
represents the forward rate , and diagram (b) the backward rate k\ ^ . Each graph has 
fixed charge order along the forward/backward path. Enqueuing the double path (a) into a 
single path yields two contributions, + © and © + , while diagram (b) yields — Q and 
— . These are taken into account by the integration limits — i oo and i oo in 163B . 



which is in agreement with the expression (|60|) . 

Consider next the conjecture (|55|) with (|4T|) and (|42|) . This yields the ex- 
pression 

Gi = -x 2 [H!(K,-p) - Hi(K, p) } . (67) 
P 

The function H± (K, p) coincides with the cumulant C 2 (p) . It is instructive that 
Hi(K,p) can be written as the imaginary-time noise correlation integral 

i fir & 2ipv 

H^K^p) = - / dv . . (68) 

7r J (2smu) 2A 

The analytically continued integrand of (|65|) . e _2pz /(sinhz) 2K ", is free of singu- 
larities in the half-strip < Rez < oo , > Imz > — tt of the complex plane 
z = u — iv. Accordingly, the closed contour integral along the edges of the 
half-strip vanishes. This entails that the noise integrals Hi(K,p) and Ji(K,p) 
satisfy the relation 

Hx{K,p) = 2e ip7r sin[(p + K)tt] Ji(K,p) . (69) 

For integer p, this relation reduces to 

Hi(K,p) = 2sm(K-K)Ji(K,p) , p e integers . (70) 



Hence the expressions (|66j) and (|67| coincide for integer p. Evidently, the result 
(|66|) of the real-time calculation is also valid for complex p. This confirms that 
the expression (|67[) , which is based on the thermodynamic method and has been 
calculated for integer p, can be analytically continued to complex p. 

Unfortunately, the real-time approach is not practical for perturbative con- 
tributions G n (x, e/T) with n > 2, since the relevant Coulomb gas multiple inte- 
grals for charges distributed on the forward/backward paths can not be carried 
out in analytic form. The loophole now is to take the imaginary time-route to 
calculate the free energy for integer p, as sketched in Section[3J With the conjec- 
ture that the resulting expression (|4Tj) with (|47|) can be analytically continued 
to complex p, we then get the conductance (|55[) . Let us now see whether this 
line of argument is correct in all orders of x. 
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4.2. The case K < 1 

In the regime K <C 1, the leading contribution to the coefficient l2n(p) in 
the multiple series (I47|) is the summand labelled with mi = ni2 = ■ ■ ■ = m n = 0. 
Furthermore, for K <C 1 the product term e,,(0) given in (|4"5)l reduces to the 
form ej(0) = (-ftT+p) ]. With these truncations the series coefficient l2n(p) 

takes the concise form 

Interestingly, the resulting perturbative series of the partition function (|4ip can 
be summed in closed form, 



r 2n+p/K 

Z TB (x,p) = T{l+p/K)x 



-p/K \" 



n — 



n\T(l + n+p/K) (72) 

= r(i +p/K)x-P/ K I p/K (2x) . 

The function I v (z) is a modified Bcssel function with index v. Applying func- 
tional relations of the Bessel functions 16[ the tunneling susceptibility dT/dx 
takes the analytic form 

—Ftb{x,p) = —T— 7 —y . (73) 

dx I p / K (2x) 



With this the expression (f5"5j) for the nonlinear conductance takes the concise 
analytic form 

Q (x,e/T ) = — 2x1m- , q = —— . (74 

9 I-iq/K{2x) 2ttT 

Upon employing recursion relations of the Bessel functions, this expression can 
be transformed into 



smh(nq/ K) 1 

iq/K 

This form reduces to the linear conductance 



g(x,e T) = 1 — — -—-j . (75) 



Gnn(x) = g(x,0) = 1- jj^. (76) 

With the substitution 2x — > Ej/T, e — >• 2eV^, and if — > p, the expression 
([71)1 corresponds to the result of Ivanchenko and Zil'berman for the conductance 
of a classic overdamped Josephson junction [l?}. For the standard derivation of 
the expression (|74p within a real-time Coulomb gas approach, resulting in an 
infinite continued fraction expression, we refer to Refs. (ill Il9l. FJ] . 

In the weak-tunneling model (|2|), the regime K <C 1 corresponds to the 
classical regime of the Brownian particle model ([I]), as the damping parameter 
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A'wb = 1/A is very large. The respective conductance is independent of A' WB 
and is 

nt irr\ i 2x t h-i q (2x) 

g MT) - 

sinh(^) 1 1 J 

as follows from the selfduality In the strict classical limit we have 2x = Uq/T 

Weak adiabatic quantum fluctuations may be taken into account (i) in the 
WB model by an effective corrugation strength Uq —> Uq = Uq (2-kT j c ) x I k . 
and (ii) in the model (fT3| by an effective Josephson energy £j 4 £J = 
E' 3 (2tvT/uj r )p 0, respectively. 

Consider next the zero temperature limit of the expression (|77p. As T — > 0, 
both the argument and the modulus of the Bessel function become very large. 
In this case the uniform asymptotic expansion of the modified Bessel function 
applies [l6| . The resulting expression for the nonlinear conductance of the weak- 
corrugation model ([H is 



G\t = o = e(e-27rtfoh/l-( — ) • (78) 



The expressions (|77|) and (|78p coincide with the solutions of the classical 
Smoluchowski equation in corresponding regimes poj . For e < 2ttUq, the con- 
ductance (|78[) is zero because the sliding motion of the overdamped particle 
comes to rest at locations with zero slope. The particle rests there forever, 
since thermal and quantum fluctuations are absent in the classical regime at 
T = 0. For e > 2ttU~o, the potential in ((T|) is sloping down everywhere, so that 
the particle is continuously sliding. 

Transcribing (|78|) into the regime of classical phase diffusion in the Josephson 
junction, which corresponds to the regime p « 1, yields for the nonlinear 
conductance at T = 



e(^-i?/ c ) 



i 



yJl-(RI c /VY 



(79) 



where I c = 2eEj is the critical current. 

In conclusion, our analysis shows that the above imaginary-time approach 
to noncquilibrium transport in the limit K — > (and K — > oo in the dual 
model) reproduces all known expressions calculated from an elaborate real-time 
method. 



4.3. The limit T -> 

The limit T — > constitutes a crucial proof whether the analytic continuation 
from integer p to complex p in the expression pij) with (|47p is correct for general 
K in all orders of x. 
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For general p, the m-fold ordered sum in (p) can not be done in analytic 
form. The simplest term (p) is a generalized hypergeometric series of the 
form 2nF 2 „-i(- ■■ ; • • ■ ; 1). 

As \p\ — > oo, the ordered sums in (p) turn into ordered integrals which 
can be solved in analytic form. With the substitution k — > pu, we have the 
mapping 



p 



du 



(80) 



and the asymptotic expansion of the (pu) emerges as 



eApu) = 



P 

T 2 (K) 

1 



<,U) 



1 



&2l(j)- 



c/ 2 



2m + Mj) 



p uU 
1 

'uU 



(81) 



where £/ = 1 + u. The expansion coefficients are 
h(j) = \K{K-l){2j-l), 
Mi) - ^A(A-l)(A-2)[3A(2j-l) 2 
622 (j) = ±A 2 (A-l) 2 (2j-l) 2 . 



1] 



(82) 



From this wc find that each term of C^ip) behaves asymptotically as 
K-2)n+ m . However, there are formidable cancellations by adding up the 
individual terms, as we immediately see from the expressions (|52[) - (|54[) . As a 
result, as far as the cumulant C2 n (p) is concerned, subleading terms in the curly 
bracket of ([BTj) become the leading ones. The analysis shows that there is an 
extra factor 

(K/pT 1 - 1 , (83) 

which results from the relevant subleading terms. Thus, in reality, the asymp- 
totic power law is '(p) oc p (2K-2)n+i^ an( j nence C 2n {p) °c p ( 2K - 2 ) n + 1 , 
Consider now first the partial cumulant coefficient 



P 



11 

C 2rAP) = -f2^jq £ 



T(K + j)T{K+p + j) 

r(i + j)r(i+p + j) 



In the asymptotic regime \p\ — > 00 we have 



P 



2(K-l)n+l r°° 



nT 2n {K) 

p 2(K-l)n+l r(1 _ n 



dit[u(l + u)] (K ~ 1)n 

nK ) T(2n - 1 - 2nK) 



(84) 



(85) 



r 2n (K) 



r(n - nK) 



The integral in the first line of (|85|) is convergent in the regime 1 — — < K < 
1— i. On the other hand, the total cumulant C2n(p) is regular down to A" = 0, 
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as follows from the integrals in ((32), i.e., < K < 1 — In comparison with 
Cznip): the expression (|85|) has n— 1 additional singularities located at if = —, 
where m = 1, 2, • ■ ■ , n — 1. The spurious singularities are contained in the term 
r(l — n+ nK). The analysis now gives that the sum of all partial cumulants in 
(|5Tj) gives C^ip) times a polynomial a n {K) which is of order n — 1 in AT, 

C 2n (p) = <*„(#) . (86) 

The polynomial a ra ( K) provides zeros just at the locations of the spurious sin- 
gularities, so that C2n{p) is indeed regular in the range < K < 1- i. With 
the limiting value a„(0) = 1, which follows from we thus get the unique 
expression 

[U ra = l (87) 

(-1)"- 1 T(nK) 



(n- 1)! T(l -n + n/\) ' 

Insertion of (|5T)) into (|55j) yields the asymptotic form of the cumulant coef- 
ficients in the regime \p\ — > oo, 

C M (-1)'- 1 T(nK) T(2n 1 - 2nK) 2(K _ 1)n+1 



In order to corroborate these considerations and the result (|88|) . we calculate 
Ci{p) and Ce(p) in |Appcndix A| explicitly. 

The cumulant C2n(p) has simple poles at K = 1 — + 2£ where m = 
0, 1, 2, • • • . These poles are due to the Gamma function term c n {K) = T(2n — 
1 - 2n/<0/r(n - nK) in C 2 „(p). 

Consider next the analytic continuation provided in the expression (|58p . 
Most interestingly, this produces a trigonometric factor cos[ — K)n], which 
has zeros just at the poles of c n (K). As a result, the perturbative conductance 
contribution G n (x, c/T) turns out to be regular for all K > 0. It may be written 
as 

_ (-i)«-i r(|)r(i + niO fq ^(K-i)n 2n 



n\ T 2n (K)T(^ —n + nK) V2 



)z 1 1\ — i ) n 
x 2n . (89) 



A crucial point now is that temperature cancels out in the combined expression 

(? 2(X-l)n :r 2n 

Finally, it is useful to combine the tunneling coupling A and cutoff u c to the 
universal Kondo scale e as [l| 



7T 2 A 2 

2-2K r,2-2K n ^» 



00) 



r 2 (iO w 2 ^ 



With use of this scale in the expression ([89)1 . the weak-tunneling series of the 
conductance at zero temperature and general K takes the concise form 

0(e,e o ,tf) - ^ —-j— (_j . (91) 

71 — 1 V ^ 7 
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This is exactly the expression for the conductance at T = found with the ther- 
modynamic Bethe ansatz by Fendley, Ludwig and Saleur [2l[ . The series (f9Tj) is 
appropriate for K < 1 in the regime e > e cr . where e cr = — K \K K '^ 1 ~ K '^ e . 

For if < 1 and e < e cr , the appropriate starting point would be the relation 
(|56]l. In this case, the resulting expression for the conductance is the strong- 
tunneling series 



g(e,e ,K) = 1- V ~ v ' y " v " ' (-) • (92) 

V ' °' ' ^ n! r | - n + n K) \eJ v 7 

Expressed in terms of the parameters of the WB model, the Kondo scale e is 



(_!)«-! r(|)r(l +n/K) /e \2(i-i/^) ra 



2-2/K _ coin 2 - 2 /* - 

Here the extra factor K 2 ~ 2 I K compared to (|90|) accounts for the mapping e — > 
e/A' in the self-duality relation (|9]). 

For A > 1 and fixed e, the regimes of weak tunneling, Eq. (|9"Tj) . and strong 
tunneling, Eq. (1921) . are exchanged. 

We conclude this subsection with the remark that there are contour integral 
representations of the conductance of which (j9~Tj) is the Taylor series and (|92| 
the asymptotic series, or vice versa 22, 23l 



5. Rates, current and noise 



The conductance Q and the current (I) = eQ/(2irK) are first moments of 
the transition rates, 



U+) _ U-) 



(94) 



Higher moments of the current give information on the statistical fluctuations 
of the transport process, 



oo 

(/M) c = n ™ + (- l ) mfc i" ) ) 



(95) 



A knowlcgde of all the connected moments (|95|) allows then a reconstruction of 
the population probability distribution. 

The dynamics of the population probability P n (t) of site n is governed by 
the master equation 



The characteristic function P(X,t) = 'Yl in & lXn P n {t) with initial state P n {0) 
S n is found from the master equation to read 



P{\,t) = ' I exp[t(e iA " - l)fcW + t (e- 



i An 



n=l 



(96) 
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The moments of the probability distribution follow from the characteristic func- 
tion by differentiation, 



( n 



">(*)) = J^n m P n (t) = (-iA) P(A,t) A = (j . (97) 



The cumulant expansion 



n — 1 



leads us to the connected moments (|95|) of the current. 

It is evident from (|95[) that knowledge of all moments is tantamount to 
knowledge of all rates fc„ in all orders of a; or A. 

The rate k^ is built up by paths which start out at an arbitrary initial 
state of the reduced density matrix (RDM), say (0,0), then may run through 
any number off-diagonal and diagonal states of the RDM, and finally end in 
diagonal state (in, in). The minimal number of moves, each of weight A/2, is 
2n. Since the paths may make arbitrary detours, the transition rate kn is a 
series of partial rates k^J , where the second index indicates that the transition 
consists of 21 moves, or 2£ charges in the charge representation, 

*4 ±} - xx? = fx****- (") 

In the last form, the dependence on the effective fugacity x is exposed. For 
paths, which have interim visits of diagonal states, the reducible part must be 
subtracted ( see Ref. 

The partial rates h^) may be calculated via analytical continuation of the 
£th cumulant of the partition function Z. In the analytical continuation the 
2£ charges representing the cumulant C2e may be arbitrarily partitioned on the 
forward and backward path, as sketched for the case t = 2 in Fig. [2] The charge 
correlation factor for m time-ordered charges {uj = ±1} on the forward path 
q(t) and t time-ordered charges {vj = ±1} on the backward path q'(t') is 

!m j-i 
UjUiQ(tj - U) 
3 = 2 i = 1 

j = 2 i = 1 



i = 1 j = 1 



x ex P { 1 6 I Y, u i *J ~ X! Ui *< 




(100) 
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where Q(i) = 2K ln[(o; c /7rT) sinh(7rTt)] + inK 8gn(t), and where the last line 
represents the correlations due to the bias. For a charge sequence contributing 

to , wc have the constraints J2j u j = ^ zn an d Yli v i = ^ zn - The term 
represents the sum of all (21 — l)-fold Coulomb integrals with the correlator 
(|100p which are compatible with the constraints. 
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+ + 
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(b.l) 

© © © - (b - 2) 
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(c.2) 

(d.l) 
(d.2) 

(e.l) 
(e.2) 

(f.l) 
(f.2) 
(f.3) 



Figure 2: The six imaginary-time charge diagrams representing Z2 or the connected part T2 
(left), and the various possibilities to partition the charges on the forward/ backward path. 
A charge on the backward path changes sign when the time axis is reversed. The encircled 
charges are the charges {v^ in J 1 00 ft . Six diagrams each contribute to and to fc[ 

and one diagram each represents feo an d ^2 ■ ^he °th er sixteen charge sequences, which 
all contribute to fcg \ are omitted. Upon enqueuing the various charge sequences along the 
forward/backward path into a single path, the graphs (a. 2) and (f.2) each split up into 6 
different time-ordered charge sequences, and each of the residual 12 graphs splits up into 4 
different time-ordered charge sequences. 
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The property Q(t — i /T) = Q*(t) in the correlator (|100p ensures detailed 
balance between the forward and backward partial rates k^J [24j |. 

*y = 6—^*45 . (ioi) 

Actually, detailed balance already holds pairwise between particular subsets 
with fixed time-order along the Keldysh contour. The subset pairs are related 
by charge conjugation and time-reversal. For instance, in order A 4 there are 
detailed balance relations between (f.2) and (a. 2), between (c.l) and (a.l), be- 
tween (e.l) and (b.l), and between (f.l) and (d.l). (cf. Fig. [3]and Ref. [24[ for 
details) . 

Unfortunately, the Coulomb integrals determining the various contributions 
to the partial rate k^J can not be evaluated in analytic form for general T and 
general K. 

However at T = 0, one may use besides the detailed balance relations for the 
subset pairs also scaling properties of the integrals because of the logarithmic 
interactions at T = 0. These together yield a multitude of relations between 
different Coulomb integrals of same order £ [24j|. As a result of the analysis, 
one finds formidable cancellations in the sum of the various Coulomb integrals 
of given order £, and only those give rate contributions at T = 0, where the 
charges {uj} and {vi} are all positive, i.e., only direct paths without detours, 

= 4|2 z 2 ", *4 _) = 0. (102) 

Since fci + ' oc x 2n at T = , it is straightforward to deduce from the pertur- 
bative series (IMj) the analytic expression 

n\ r(f -n + nK)2ir \ eJ y ' 

Expression (|96|) with (|103j) yields an analytic form of the characteristic function. 
In addition, we get the concise moment relation 

<'<™'>c= (5 S) CO- P04) 

Thus, all statistical properties of the quantum transport process at T = can 
be deduced directly from the current [25|. 

Let us finally return to finite temperature. Evidently, the analytic continu- 
ation rule in Eq. (|55[) can not provide information about individual tunneling 
rates. Rather we find with use of (|M|) and 

00 £ 

Im [F{x,p) - F{x, -P) ] = £ £ 7 (fcS " *#) • (105) 

l=ln=l 

This yields for the term of order x 2t 

Tx™ Im [C 2l (-V) CM 1 = £ 7 (^S - (106) 

n = l 
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Thus, the analytically continued perturbative series of the free energy yields for 
each term, say Tn, a particular linear combination of partial rates of same order 
x 21 , as given in Eq. (|106[) . In order to calculate individual partial tunneling 
rates, which are required in the full counting statistics at finite temperature, an 
analysis of the Coulomb integrals in real-time is indispensable. 

6. Conclusions 

Wc have discussed a variety of seemingly different physical models which 
display a wide range of interesting characteristics and can be realized exper- 
imentally. These models have in common that they can be mapped on a 
one-dimensional quantum field with gapless bulk excitations and a sinusoidal 
boundary interaction. The masslcss boundary sine-Gordon model represents 
an impurity in a Luttinger liquid. When the Luttinger liquid modes are in- 
tegrated out, the model describes a quantum-dissipative particle moving in a 
tilted sinusoidal potential. All these models can be successfully treated by a 
variety of interesting and powerful techniques. In most cases, the task of cal- 
culating nonlinear transport and full counting statistics requires application of 
noncquilibrium techniques. For general model parameters, this is a challenging 
endevour. Despite the fact that the systems are integrable with help of the 
thermodynamic Bcthc ansatz, this method can not be used to calculate the full 
counting statistics, apart from the current. To gather full information about 
the statistical fluctuations, a real time calculation of the individual tunneling 
rates is essential. 

Here we pursued a different route. We seized a suggestion by Fendley, Lud- 
wig and Saleur, who calculated the twisted partition function of a "log-sine" 
Coulomb gas on a ring, and analyzed a conjectured relation with the nonlin- 
ear conductance which involves analytic continuation in the bias and coupling 
strength. We studied the weak-tunneling limit, the classical regime, and the 
zero temperature case and found in all greement with results available 

from a non-equilibrium real-time calculation. 

Appendix A. The cumulants Ci(p) and Ce(p) 

Appendix A.l. Ci(p) 

With the substitution ([50]) and use of the asymptotic expansion (|5Tj) . the 
expressions (|53[) and turn into 




(A.l) 



(A.2) 



where 



/(«) 



it 



K-l 



(1 + U) 



K-l 



(A.3) 




(A.4) 
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The function h(u) originates from the terms of order 1/p in the curly bracket 
of (f8Tj) . The integrals (|A.1[) and (|A.2[) are convergent in the regime \ < K < |. 
Observing that 

f(u)h(u) = Kf'{u) , (A.5) 

where the prime indicates differentiation with respect to the variable, the double 
integral can be reduced by partial integration to a single integral, which is a Beta 
function. In the end, we get 

cMftO - i P ' K - 3 r(2K i)r(3 - 4g) 

° 4 [P) ~ 2 T*(K) T(2~2K) ' (A ' bj 

c| 2) (p) = -2^C| 1) ( P ). (A.7) 

This yields a 2 {K) = 1 — 2K. Thus, in the sum of C±(p) and C^\p) the singu- 
larity at K = 1/2 is abrogated, and we end up with the cumulant expression 

U[p) ~ 2 T\K) r(2 -2K) P ' (A - 8) 
which is regular in the range < K < |. 

Appendix A. 2. Cs(p) 

As \p\ — > oo, way may use ([80)) and the asymptotic series (|8Tj) in the expres- 
sions (|54|) . We then obtain 

C& {p) = 3fwi) duf{uh (A ' 9) 

6 A'— 5 r roo nu 

C 6 2) (P) = 3^)[y o du(f(u))'J o dvf(v) (A.10) 

poo p u 

+ / d«/'(«) / du/ 2 H 

Jo Jo 

, \ n 6K-5 



duf"{u) / dvf(v) / dw/M 
Jo Jo 

+ 2 I duf'(u) / dvf'(v) / dwf(w) 



Here we have again employed the relation (|A.5[) . The integrals in (|A.9|) - (|A. 1 1|) 
are convergent in the region | < K < |. The double integrals (|A. 10|) and triple 
integrals (|A.11|) may be transformed by partial integrations into single integrals, 
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which again are representations of the Beta function. We get 



(i)^ _ 1 P 6K ~ 5 T(3K - 2)r(5 - 6K) 
3T 6 (K) T(3-3K) 



= ~K(f>(p), (A.12) 



The sum of these terms yields Cq(jj) = (Xz(K) Cg (p) with 



a 3 (K) = ±(3K-l)(3K-2). (A.13) 

~\ nr , n + „A »t lc' — 2l 



The prefactor (|A. 13|) just cancels the singularities of C 6 (p) located at K 



and at K = i. While C@ (p) is regular in the range | < K < |, the resulting 
cumulant expression 

C W -*(IOCfW-l^^^-, (A.!4) 



is regular in the extended regime < K < 
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